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ABSTRACT 


Short and long term variations in the solar constant are examined theo- 
retically. The variations observed by the Solar Maximum Mission, lasting 
several days and associated with the passage of sunspot groups, strikingly 
demonstrates the well known lack of a "bright ring" effect around sunspots. 
This suggests that sunspot magnetic fields do not simply block the heat 
flowing upward into the photosphere. Rather, it is suggested that gravita- 

tional draining occurs; this cools sunspots and transports downward the heat 
that would otherwise flow into the photosphere. A model of sunspot 
temperature with depth shows modest support when compared with the empirical 
model of Van't Veer. Secular trends in the solar constant may occur and be 
associated with the influence of the convection zone magnetic field upon 
convective heat transport. As a start to understanding this problem, the 
Schwarzschlld criterion has been modified to include the effects of magnetic 
field. 

INTRODUCTION 


Recently investigations by Livingston (ref. 1), Kosters and Murcray (ref. 
2), Livingston et al. (ref. 3), Dicke (ref. 4), and Willson et al. (ref. 5) 
suggest that short term and secular changes in solar luminosity may be 
occurring. Theoretically, the influence of the solar cycle magnetic field on 
solar luminosity is a multifaceted question, since there are numerous ways in 
which magnetic field can affect the sun's luminosity. 

Thomas (ref. 6) has theoretically investigated luminosity changes 
associated with effective solar radius changes produced by magnetic buoyancy 
in spots. Hoyt (ref. 7) looked at differences in luminosity associated with 
observable photospheric features-spots, faculae, etc. - to obtain a measure of 
possible global luminosity differences. He found a relationship between 
umbra-penumbra area ratios and terrestrial temperature variations. Willson et 
al. (ref. 5) have examined dips in the solar luminosity associated with the 
passage of spot groups past central meridian. These dips are a "short term" 
influence %ie shall discuss shortly. "Long term" influences may arise from 
changes in the convection zone. Schatten (ref. 3) suggested that magnetic 
buoyancy may influence convection; Splegal and Weiss (ref. 9) argue that the 
heat transport by convection can be affected by the magnetic field, and that 
the influence is most significant at the base of the convection zone. We 
shall discuss, later, a mechanism that allows the- long term influence of the 
magnetic field upon the convection zone to be calculated from stellar 
models. First, let us examine the short term variations associated with spot 
groups. 
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SHORT TERM VARIATIONS ASSOCIATED WITH SPOT GROUPS 


As the photospheric material cannot serve as a reservoir for the vast 
flow of energy on time scales longer than about a second, a "bright ring" 
around sunspots might be expected if the Biermann (ref. 10) mechanism 
involving field inhibition of convective heat transport occurred. 

A new school of thought developed, beginning with Danielson (ref. 11), in 
which the convective generation of Alfven waves below the sunspot consumed so 
much energy that cooling of the spot occurred with an attendant field concen- 
tration. It also eliminated the "bright ring" effect wherein the blocked heat 
would flow around the spot into the surrounding photosphere. The lack of a 
bright ring was strikingly demonstrated by the recent "solar constant" obser- 
vations of Willson et al. where dips in the solar constant occurred when 
sunspots approached central meridian. If any bright ring were present, it 
would cancel the spot energy deficit and no dip should occur. It is these 
dips that are the subject of the present section. Other Investigators have 
also found changes in the solar constant with sunspot visibility. Foukal and 
Vernazza (ref. 12) found a level of 3 x 10"^ change in the solar constant 
(similar to the SMM findings) and Chapman (ref. 13) found a dip in the sun's 
brightness equal to 62% of a sunspot's area. 

Here we have followed much of the theoretical guidance suggested by Meyer 
et al. (ref. 14) and by Parker (ref. 15, 16). We further this third view that 
pores, knots and sunspots are a kind of dynamical solar sink in which material 
drains gravitationally downward. This utilizes a dynamical gravitational 
draining as the mechanism to explain 1) the cooling of features (ref. 15) 
2) the elimination of the bright ring around sunspots, and 3) the field 
concentration mechanism (ref. 14). It enables a calculation of the temper- 
ature of these downdrafts to be roughly 1000°K cooler than the photosphere 
from basic principles and that the downflow in fluxtubes should be roughly 
several km/s. 

Parker (ref. 16) has suggested that the "6-8 km s~^... downdrafts in the 
fluxtubes (leading to) concentration to 1500 G may be a direct dynamical 
consequence". This downdraft inferred by Deubner (ref. 17) from observations 
may also be Important in stabilizing fluxtubes Insofar as the magnetic 
pressure of the fields, in a solar atmosphere with no inward and downward 
flow, would tend to disperse the fields in a short amount of time. Meyer et 
al. (ref. 14) discuss the stabilizing effect of this inflow of material as a 
hydrodynamic "collar" which the spot wears. 

We discuss here, that this inward and downward flow as shown in Figure 1, 
may also play a role in the cooling mechanism for spots, as well as pores 
where Frazier (ref. 18) has observed downdrafts up to 3 km s**^. For pores and 
knots, large velocities in the photosphere are found, or at least inferred 
however, for spots Reciters (ref. 19) notes that they "show no vertical motions 
... exc'iedlng 25m/sec." Some confusion has arisen "because of the limb effect 
of the surrounding photosphere, sunspots appear to have a downflow of 
~400 m/sec ...." Thus do%mflows are seen or Inferred for small flux tubes, 
and not for spots where we hypothesize a deeper gravitational draining occurs 
owing to the larger size of the magnetic object. In this view, the field 


208 



originates from dynamo processes, the final stage of which Is magnetic 
buoyancy, whereby a large concentrated field erupts at the sun's surface. In 
the absence of a suitable cooling and field stabilisation mechanism, the 
strong field would quickly dissipate. However, as the material near the sun's 
surface Is continually radiating luminous energy Into space, the gases cool 
and with their Increased density return In conduits to the heat source from 
which they came, thereby completing the convection cycle. Thus by analogy 
with any thermal cycle, the material i ovlng towards the heat source would be 
the coolest. In the absence of sunspot fields, the cooled photospherlc gases 
return at the boundaries of the supergranulation, and so form pores. In the 
presence of sunspot fields, they are aided in their return by this field 
conduit, which tends to reduce the turbulent viscosity and so provides an 
easier pathway. Owing to the larger size of a spot, not enough material 

appears capable of congregating Into an observable downflow in the photosphere 
and we hypothesise It does so at depth. 


To calculate the size of the effect, we take a volume element within the 
magnetic region, as shown In FlKurf* 1. It Is bounded by 4 sides (1-4) with no 
Interaction on the 5^" and 6^" sides due to azimuthal symmetry. He shall 
simplify the picture by assuming the following. The flow is roughly Inward to 
the box at 1, and outward at 4 . If the photosphere were not yet cool, region 
2 would radiate the same energy rate, F, as the remainder of the photo- 
sphere. Due to the ln\ nrd and downward flow, little or no convective energy 
is being transported into the region through 3 (where In the normal 
photosphere, a rate, F^ Ealances the outflow through the top). For the 
purpose of tbs calculations, the vulunte extends from the photosphere t >■ 1 to 
roughly t • 3 so that most of the radiant flux, F ■ 6.4 x 10^® erg cm“^s"* 
(ref. 20), Is Indeed emitted from the volume through 2. Further, the material 
flotis Inward through 1 at several km sec*^ and out through 2 at the same rate 
and pressure, so that we can Ignore compressive heating or expansive cooling 
as a consideration. Then the temperature difference can be shown to be: 


AT - 


FL 


“ 1,100 K 


( 1 ) 


5 p hvR 

where R is the gas constant, L and h are the length and height of the voluM, 
taking L/h > i and v - 7 km s~^. Thus for a spot, an umbral temperature of 
4,600^K is calculated showing good agreement with observed umbral 
temperatures. Put simply, the mechanism can provide adequate cooling for 
pores, magnetic knots and sunspot umbra. We now examine the question of the 
downflow. The cooled photospherlc gases, being denser than the hotter 
photospherlc gases, would tend to descend. Will this gas reach a terminal 
velocity of several km sec"^ as it descends? We consider a calculation for 
the terminal velocity of a blob of gas In a viscous stratified atmosphere. 


As the terminal velocity la governed by the Interaction with the 
stratified medium, we choose the scale height, H, to be a characteristic 
length dimension for consideration. The downward force upon the gas Is: 

F - g Ap (2) 
where g la the acceleration of gravity - 2.7 x 10^ cm s~^, Ap Is the 
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AT 1 lOO'^K 7 

additional density of the blob of gas • “y P — 4 x 10"' gn cm’-’ - 0.7 

X 10”^ gm cm"^, and H is the scale height * The viscous balancing 

upward force is: 

F - a nvH (3) 

where n is dynamic viscosity, pv, and v is the kinematic viscosity 1/10 vl. 
Busse (ref. 21) obtained a value of 10*^ cm^ s"^ for v, yielding 3 x 10^ gm 
cm'U'"^ for n. The quantity v is the velocity difference between the blob and 
the surrounding material, and H is the scale height of the photosphere 
again. Equations 2 and 3 can be solved for v, yielding: 

V - gApH^/n. (4) 

This gives 7 km s”^ for the terminal downward velocity of gases, a sufficient 
downflow to provide cooling. It is comparable with the 2 km s“* velocity that 
Parker (ref, 15) required for a downflow tdthin spots to account for their 
cooling. 

The lack of a bright ring around the gas central to the problem of solar 
constant dips is understood in this model by examining Figure 1, As the 
upward heat flux, F, transported principally by convection, encounters the 
sunspot magnetic field, the energy and gas are entrained amongst the 
descending gases where the energy is carried down toward the base of the 
convection zone and so is lost to the sunspot umbra. 

Table 1 shows the temperature and density of the umbra and surrounding 
photospherlc material. The photospheric parameters are shown from the solar 
model of Endal and Sofia (vef, 22), The umbral parameters were obtained from 
a computer model having been calculated utilizing the simple theoretical model 
outlined. That is, there is no convective heat transport into the spot from 
below and the photospherlc material is cooled by radiation into space. 

For comparison, the semi-empirical model of Van't Veer as outlined in 
Tandberg-Hanssen (ref, 23) is shown in table 2. It is based on a method 
developed by Van't Veer that determines the umbral parameters using the 
measured intensity In the wings of certain Fraunhofer lines. Near x - 4, the 
values agree remarkably wellj however, near x ■ 1 the agreement is less 
good, 

LONG TERM VARIATIONS IN THE SOLAR CONSTANT 


We consider, in this section, the effects of the solar cycle magnetic 
field embedded within the convection zone upon the solar luminosity. We 
develop a version of th.s Schwarzschlld criterion which includes the effects of 
magnetic field. We consider the approximation that the magnetic field 
provides a net isotropic pressure. 


The Schwarzschlld criterion (ref, 24) is obtained by considering density 
changes associated with a rising convective element of the star, leading to 
Instability If the radiative temperature gradient exceeds the adiabatic 
gradient: 
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If we now Include a magnetic field, there will be a magnetic pressure Pg, 
associated with a rising bubble. In any discussion of the magnetic stresses, 
we may neglect the tension term In the stress tensor because the field Is 
continuous and divergence B eqtials zero. Thus across any small element, the 
field tension balances on either side and Imparts no force to the gas. The 
total pressure P^, equals the gas pressure plus magnetic pressure: 

where 6 Is the gas to magnetic pressure ratio. The adiabatic density gradient 
now becomes; 

d In _ d In P^ _ d ln(H^~^) ^ 
dr •'a " y dr ■'a Y dr ■'a 

where we have used the adiabatic relation, P.V " constant, 
gradient Is similarly modified: 

d In p- > d In P^ .^ _ d In T -> _ d Ind-t^g"^) ^ 

dr ■'r " dr ^•'r ” dr ■'r ” dr -*r 


(7) 

The radiative 
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Instability occurs If the adiabatic density gradient has a larger 
magnitude than the radiative gradient. Frm equations 7 and 8 this Implies: 

A In T v-1 ^ ) dln(l+B ;^) 

lilLll > III + 2 S— f9^ 

d In P^r ' Y y d In P_ d In P_ ^ 


where adiabatic and radiative subscripts on 8 allow for differing values 
between convective bubbles and the surrounding environment. Equation 9, our 
modified Schwarzschlld criterion, can be rewritten as: 

. t- , . . , d ln(l+B:b 
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where we have Incorporated the two expressions on the right of 9 into one 
through the utilization of a parameter, 5. Here 5 equals one under the assumed 
approximation, for the onset of Instability. For the more general case 6 
depends on the ratio of field Inside the rising gas bubble to field outside, 
and geometry. Other formulae can be developed to Include the effects of other 
geometries (see Thomas and Nye (ref. 25 1. Further, If one excludes the 
Interaction of the bubble with the material, an adiabatic equation of state 
can be written, but we have allowed this to be Incorporated Into 6. As we have 
developed our criterion for the onset of Instability when the field Is not 
perturbed, we have Ignored the possible restoring force (or tension) between 
the base field and a rising bubble %ihlch may ensue with turbulence. The above 
criterion assumes that once the magnetic field and temperature gradient are 
suitable for Instability, the Raylelgh-Taylor Instability or another plasma 
Interchange Instability %flll develop to allow the bubble to rise. 


The term on the right of equation 10 modifies the adiabatic condition In 
the following fashion. Let us consider 8“^ to approach 1 In a region of the 
sun and examine how this affects the onset of Instability. Figure 2 shows the 
magnetic field, the ouignetlc to gas pressure ratio (8~^), *nd the term in 
equation 10 modifying the Schwarzschlld criterion. As the gas pressure 
decreases radially outward, the whole term In the square bracket will first 
exceed one, and then fall below one, to possibly a negative value. Thus a 
single layer or region of magnetic field will form a relatively stable layer 
In Its lower side and a relatively unstable layer In Its upper side. The word 
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relative ia used to indicate that these layers may be stable or unstable, 
depending upon the remainder of the Schwarsschild condition, and thus, they 
are merely isodifylng the condition. However, in case the modifying term 
reaches a negative value, as the left hand side of equation 10 is positive 
definite, instability is required. This is conventionally referred to as 
magnetic buoyancy, and in the usual treatment only the field considerations 
are provided. We see here, hotiever, that ordinary convection may also be 
affected by the presence of a subsurface field. 

If we consider the fields in the solar Interior suggested by dynaiao 
theory (ref. 26), the term on the right of equation 10 appears to be sig- 
nificant only at the - base of the convection zone (except in centers of 
activity). If the gas pressure variation has a scale height Hq, and the 
magnetic pressure a scale height Hg, the term on the right can reach a 
magnitude of order Hg/HQ* One of the problems associated with solar dynamo 
theory is the formation of a region where magnetic fields may regenerate 
without being lost to magnetic buoyancy (ref. 26, 27). In this result the 
underside of a magnetic field region tiould be a stable location where fields 
could regenerate. Above the region, if the field approached a high enough 
value, instabilities could form with a balance developing between regeneration 
and deterioration. Unno and Ribes (ref. 28) have examined magnetic buoyancy 
and have found that 200-300 gauss dynamo field can be retained in the con- 
vection zone for time scales up to about 80 years due to turbulent viscosity. 
The large scale field, although retained, could affect the onset of convective 
instability through the above criterion, thus the process of magnetic stabil- 
ization and destabilization may influence heat transfer at the base of the 
convection zone with solar cycle phase much as Spiegel and Weiss (ref. 9) 
suggest. 

CONCLUSIONS 


Downflow of material within sunspots is seen as a key to the 
understanding of their stability, temperature and heat flow. The Inward 
velocity prevents the magnetic field from expanding Into the surrounding 
surface hy balancing the spot's magnetic pressure against the velocity change 
from an Inflow into a downflow. The cool, dark characteristics of sunspots 
are to be explained by their being the cool end of a thermal cycle. In this 
cycle, the gases have radiated their energy Into space, are allowed to descend 
along field conduits into a cooler, denser state and cut off the heat flow. 
The heat flux that would normally be convected into and radiated by the 
sunspot Is carried downward by the flow to deep layers In the convection 
zone. Thus, the minima in the solar output associated with sunspot groups 
observed by Willson et al. (ref. 5) are explicable. 

Further, for secular changes, we consider the magnetic field at the base 
of the convection zone to affect the Schwarzschlld criterion and form 
stabilising and destabilizing layers that may cause a significant variation in 
the solar luminosity similar to that in the model of Spiegel and Weiss. 


& 
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Table 1 


Depth 

Density 


T 

^ u 

km 

10“^gm cm”^ 

lO^K 

lO^K 

0 

1.69 

5.75 

4.60 

kl 

1.94 

6.46 

5.04 

98 

2.08 

« 

00 

5.28 

238 

2.98 

10.13 

6.94 

448 

5.18 

11.57 

9.73 


Table 2 


t(5000) 

Depth, km 
from t*.5 

lO^K 

T 

u 

10^°K 

.01 

-220 

4.70 

3.57 

.05 

-135 

4.93 

3.70 

.1 

-100 

5.07 

3.78 

.3 

-50 

5.51 

4.02 

.5 

0 

5.83 

4.18 

1.0 

40 

6.41 

4.48 

2.0 

60 

7.18 

4.84 

4.0 

100 

8.10 

5.24 
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CROSS SECTION OF A SUNSPOT 


Figure 1. Shown is the field geometry and downflow velocities 
in the vicinity of a sunspot or pore. The inward and down- 
ward flow beneath the photosphere and below acts as a "collar," 
stabilizing the sunspot. This flow also returns the cooled 
gases to the base of the convection zone. It is by the liber- 
ation of radiant heat that the gases are cooled to umbral 
temperatures. The energy flux that would flow upward into 
the sunspot is entrained within the cooler gases and descends 
to the convection zone base to be reheated. The box at the 
top labelled with sides, 1-4, exhibits a simple heat calcu- 
lation for tt sunspot. 
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Figure 2. Schematic showing the radial iation of magnetic 

field B, magnetic to gas pressure ratio and the term 

modifying the Schwarzschild criterion. As can be seen when 
the modifying term exceeds 1, the instabil'cy is more diffi- 
cult to attain as convective elements have difficulty rising 
through a more tenuous field saturated layer. On the upper 
side of a field layer, however, destabilization occurs, and 
convective elements may more easily rise. When the modifying 
term becomes negative, then absolute buoyancy occurs as the 
field strength is sufficient to cause instability even 
against an unfavorable temperature gradient. 
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